Influence of statistical fluctuations on K/n ratios in relativistic heavy ion collisions 
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The influence of pure statistical fluctuations on K/n ratio is investigated in an event-by-event way. 
Poisson and the modified negative binomial distributions are used as the multiplicity distributions 
since they both have statistical background. It is shown that the distributions of the ratio in these 
cases are Gaussian, and the mean and relative variance are given analytically. 
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The ultimate goal of current and future ultra-relativistic heavy ion program is the production and characterization of 
an extended volume of deconfined phase of quarks and gluons, the quark-gluon plasma (QGP). The possible existence 
of QGP as the equilibrium state of strongly interacting matter has been predicted in lattice theories jjj at sufficiently 
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^vq . high temperature and density. So QGP might exist in the early universe and inside the neutron stars. Heavy ion 
collision at extremely high energies is the only chance to study QGP in a controlled way. In laboratory experiments, 
O however, the QGP can survive only for very short period of time, and it will cool during its expanding and become 
hadronic matter at freeze-out point. It is very difficult to detect QGP in experiments because most events are, due 
to collision geometry and quantum fluctuations, without QGP even if the condition for QGP creation can be reached 
in the experiments. So, a variety of possible signatures for the transient existence of the deconfined state of matter 
in the collisions have been proposed theoretically and experimentally (see || Q for a review). One of the most 
I , important and most interesting signatures is the enhanced strangeness production in the process of heavy ion collisions 
' H compared with that in p — p collisions. The investigation of strangeness production is useful for understanding the 
mechanism of heavy quark production since no strangeness content exists in the initial state in heavy-ion collisions. 
C • Strangeness production can also shed light on the time scale of chemical freeze-out and therefore can carry information 
C\ about the early stage of a heavy ion collision [||. 

Great progress has been made with the proposal and applications of event-by-event analysis Jj| |[ ||, [n], |ll|, [j^, [l3| in 
the study of high energy heavy ion collisions. Event-by-event analysis of heavy ion collisions became possible with the 
advent of large acceptance detectors jbi] . The philosophy of event-by-event physics is based on following speculation: 
Although the conditions to produce QGP may be reached in every event, the fact that a phase transition is a critical 
phenomenon implies that it may occur only in a very small sub-sample of events. So the fluctuations accompanying 
II 1 the phase transition will, in effect, be averaged out in the conventional ensemble analyses. The event- by-event analysis 
searches for fluctuations of observables at the event level, so it can be used to select interesting or anomalous event 
candidates ]l5| with specific dynamical properties. So the new method can provide dynamical information which 
cannot be obtained from the traditional inclusive spectra. 

Recently the event-by-event analysis is used in experiments to analyze strangeness productions in high energy heavy 
ion collisions. In |l6|, [lTj the kaon to pion ratio for single event is studied as the event observable with which one 
can look at fluctuations in the chemical freeze-out stage of the collisions. Of course, the distribution of the K/n ratio 
is the result of a combination of statistical and possible dynamical fluctuations in the strangeness production. The 
method of mixed events, in which there are no momentum correlation among particles by construction, is used in the 
experimental analyses to estimate the effect of the former on the distribution of K/n ratio. Both distributions for 
K/n ratio from the real data and the mixed events are Gaussian, and the widths are only with a little difference. That 
means the dominance of statistical fluctuations in strangeness production. Although it is well established that there 
is no spatial-temporal correlation in the mixed events, we do not know exactly if there exist some other correlations. 
In the study of K/n ratio the only quantities concerned are the numbers of produced kaons and pions. The event- 
mixing technique will not change the total numbers of the particles produced. These numbers carry some important 
information about the chemical fluctuations at the freeze-out point. So, some information about the strangeness 
production may still be retained in the mixed events. To investigate the problem theoretically further one should ask: 
What is the distribution of K/n ratio if there exist only statistical multiplicity fluctuations in the collisions? Will it 
be a Gaussian as showed from the mixed events? If yes, what are the mean and the width? Are they in agreement 
with those obtained from the mixed events? What will be the dependences of them on mean multiplicities of kaons 
and pions? 

In this paper the K/n ratio for the case with pure statistical fluctuations is studied in an event-by-event way with 
two assumptions: (1) There is no dynamical fluctuation in the productions of both the kaons and pions, and (2) there 
is no correlation between the productions of the two kinds of particles. This is, of course, a trivial case, but it is the 
base for any event-by-event investigation on the influence of dynamical fluctuations on K/n ratio. 

Since we have assumed no dynamical fluctuations in the productions of kaons and pions nor correlation between 



o 
o 

Q 
i 

O 



X 



2 



them, the distributions of the multiplicities of kaons and pions can be specified from some statistical considerations. 
Two possibilities will be considered in this paper. As a first possibility, the distributions can be Poissonian, i.e. the 
probabilities with which there are k kaons and/or it pions in an event are, respectively, 

Px = fe"W, Pn = ^e-W (1) 
k\ ir\ 

with (k) and (it) the mean multiplicities of kaons and pions in the collisions. In a global analysis, the global K/ir 
ratio is then 

R g = (k)/(n) . (2) 

From the view point of the theory of coherent state the assumption of Poisson distributions for the multiplicities 
implies that both kaons and pions are emitted from some coherent sources and that no thermal effect is taken into 
account. Now we turn to event-by-event analysis of the ratio. For an event with k kaons and ir pions the K/ir ratio is 

r B = - , (3) 
and the probability density for K/ir ratio to be in neighborhood of r e can be given by 
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because the productions of kaons and pions have been assumed to be independent. 

Since px and pu are peaked at (k) and (ir) respectively, it is easy to see that P Tc will be peaked near R g . A natural 
question is: Can P re be Gaussian for given (k) and (tt) ? Since Eq. (||) is too complicated, Monte Carlo simulation 
can be used to answer this question. To get an event for this case, one needs only to generate two independent 
random integer numbers (according to two Poisson distributions with means (fc) and (tt} respectively). The random 
integer number corresponding to (fc) is taken to be the kaon multiplicity in the event, the other for pions. Then one 
gets a corresponding K/ir ratio r e for the event. The distribution of r e from the simulation is given in Fig. 1 with 
1 million events used. In the simulation the mean multiplicities for kaons and pions are chosen to be 38 and 200 
respectively for illustration. Here we choose R g = 0.19, the same as shown experimentally in Pb-Pb collisions at 158 
A GeV, for later comparison with experimental result. One can see that the distribution of the K/ir ratio is indeed 
approximately a Gaussian. This is consistent with the conclusions given from the mixed events in Refs. |l(], [l7|. To 
get more information about the distribution one may want to know the mean and width of the distribution. From 
Eq. (Q) one can see that the distribution of r e is not symmetric about R g , so the mean value of the distribution is, to 
some extent, different from R g . The mean value of the distribution for the trivial case can be calculated analytically 
as 

, oo oo , / oc oo 

R ° = (") = EE-Mfc)Pn(Tr) / £X>x(fc)l>n(*) (5) 

fc=0 7T=1 / fe=0 7T=1 



7T=1 
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In the transition from Eq. (j5|) to Eq. (||) condition (tt) 1 has been used. The relative accuracy in Eq. (^) is of 
the order of e _ ' 7r ', thus the approximation is quite good for (n) a few hundreds which is quite normal in current high 
energy heavy ion collisions. By using the generating function for the Poisson distribution 

Gn(x) = £ HLf_e-W = £ x> n (.) = e<*><- l > , (7) 

7T=0 7T=1 

the mean value of the distribution R e can be expressed as 

R e = R g j\ x{ ,) G ^- x G ^ . (8) 
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We are particularly interested in the ratio R between R e and R g , which is a measure of the deviation of the two K/tt 
ratios obtained from global and event-by-event analyses, respectively. In our trivial case 
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(9) 



which depends only on the mean multiplicity of pions. It can be seen that R will approach 1.0 when (tt) — > oo. In this 
limit the mean values of K/tt ratios from global and event-by-event analyses are equal, as can be naively expected. 
The change behavior of R as a function of (tt) is shown in Fig. 2. R is large for lower mean pion multiplicity and 
decreases with (tt). But the values of R are always close to the trivial value 1.0 within an accuracy of about 1% in 
the graphed range of the mean pion multiplicity. So one can hardly find the difference between R e and R g . Using the 
generating function for pion multiplicity distribution, one can also calculate the width of the K/tt ratio distribution. 
First one can get 
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where (k 2 ) = (k) 2 + (k), since the distribution of the kaon multiplicity is assumed to be a Poissonian. The width of 
the K/tt ratio distribution is then 
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The relative variance is then 
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Different from R the relative variance depends on both (tt) and R g (thus on (k)). The behavior of the relative variance 
as a function of (tt) for given R g = 0.19 is shown in Fig. 3. The value of R g is chosen for illustration and is the same 
as that obtained from current experiments on Pb-Pb collision at 158 A GeV. For (tt) from 100 to 300 the relative 
variance decreases quite quickly from about 28% to about 16.5%. The relative variance given experimentally from 
the mixed events is about 23% at (7r}=270.13 [jl8|. Since the relative variance depends on R g , different curve for the 
dependence of a/R e on (tt) will be obtained if another value of R g is taken as input. For fixed (tt), the larger R g , the 
smaller a/R e . The decrease tendency of a/R e with (tt) will be similar for all R g . The decrease of the relative width 
a / R e with the increase of (tt) can be anticipated from the fact that both the distributions for kaons and pions are of 
Poisson in which the relative widths of the multiplicity distributions are inversely proportional to the root of the mean 
multiplicities. In Refs. [Iq, [l?]] the mean multiplicity of positively charged pions is 270.13, in the range discussed in 
this paper. The experimental point is given in Fig. 3 by a solid square. One can see that the theoretically calculated 
relative width of the K/tt ratio distribution from the Poissonian multiplicity distributions would be too small. The 
discrepancy of calculated a j R e from the experimental result for the same (tt) shows the important role played by 
the decoherent effect in the emission of the pions and kaons. The decoherent effect will broaden the multiplicity 
distributions of kaons and pions, so that a larger relative width of the K/tt ratio distribution can be expected. Such 
thermal effect cannot be washed out in the so called mixed events. 

As the second possibility for the multiplicity distribution with pure statistical fluctuations, one can work with the 
modified negative binomial (MNB) distribution Jl9], ^0|. As a generalization of the negative binomial distribution, 
the thermal effect plays an important role in the MNB distribution. In this case, the generating function is 
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l-r(x- 1) 



(13) 



with three parameters A, r and k. As shown in the first paper in p9| , k can be interpreted as the maximum 
number of fireballs or clusters in some initial states, A and r can be related to the rates of birth and immigration 
processes. The mean multiplicity in the MNB is (n) = k(r — A). The width squared of the multiplicity distribution 
is w(n) = (n 2 ) — (n) 2 — (n)(l + r + A). With the MNB distributions for multiplicities of kaons and pions, the 
distribution of event-by-event K/tt ratio is shown to be also a Gaussian, similar to but broader than that given in 
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Fig. 1 . Our calculations show that the mean value and the width of K/tt ratio distribution have weak dependence on 
the choice of k. Thus, the parameter k is chosen to be 200 in the following. In the calculation of the dependences of 
the mean value and width of the K/tt ratio distribution on the pion mean multiplicity, there are two free parameters 
and Wk = ((k 2 ) — (k) 2 )/ (k) which represents the relative width of the distribution of kaon multiplicity. Then, R 
can be calculated from Eq. with Gmnb in place of Gn . The width can also be calculated as 
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( 1 + f 1 dx G mnb(:z)-Gmnb(0) j 1 

^UJl f *_1 , (14) 
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We do not write Wk in Eq. (|12|) as it is simply 1, because there the kaon multiplicity distribution is assumed to 
be a Poissonian. Data from NA49 experiments show that ((N 2 ) — (N) 2 )/(N) is about 2.0 for all charged particles. 
Considering the fact that most (about 90%) charged particles are pions, one can fix approximately «i, = ((ir 2 ) — 
( 7r ) 2 )/( 7r ) = 2.0 for comparison with the experiment. Numerical calculations show that a/R e has a quite weak 
dependence on w % . From the MNB distribution one can get the two parameters in Eq. ( jl3| ) for pion multiplicity 
distribution as r = (id, — 1 + (7r)/fc)/2, A = (w w — 1 — (n)/k)/2 with k preset to be 200. Dependence of a/R e on the 
global ratio R g can be expected from Eq. ( fbl] ) in the same way as from Eq. (12), i.e., the larger R g , the smaller a/R e 



for fixed {tt) and other parameters. It is shown that a/R e increases considerably with u>k- As argued in Ref. pi| , 
the main contribution to the fluctuations of the K/tt ratio comes from the fluctuation of multiplicity of kaons which 
are the particle species with fewer mean multiplicity. Normally, the relative multiplicity fluctuation for such particle 
species is larger. So in following calculations, three values of wk = 2.0, 2.5, 3.0 are chosen to present the dependence 
of a I ' R e on the width of kaon multiplicity distribution. Because the mean ratio is independent of Wk, the calculated 
R are the same and close to 1.0 within 2.5% with the three choices of wk, and are given also in Fig. 2. However, 
er J R e increases considerably with wk, as shown in Fig. 3. In those analyses in |l6|, |l7], |l8| cr/R e from the mixed 
events is about 23% with (7r)=270.13. Then from Fig. 3 one can see that the calculation from the MNB distribution 
can give value of a/R consistent with the experimental result from the mixed events when Wk is suitably chosen (less 
than 2.5). 

However, the value of wk in real experiment is unknown yet and may be different from 2.5 significantly. Thus, 
the calculated result from pure statistical consideration cannot at present be compared directly with that from the 
mixed events out of experimental data. The possible difference between them may have deep physical implications. 
This may show that there exist some other physical effects in the event-by-event K/tt distribution, which might 
be the correlation between the productions of kaons and pions. Such effects should be studied theoretically and 
experimentally. 

As a summary, the statistical influence on the K/tt ratio is studied in an event-by-event way. It is shown that 
the distribution of the ratio in the pure statistical case is Gaussian with mean and width depending on the mean 
multiplicities of kaons and pions. It is interesting to further study the effect on the ratio from correlations between 
the productions of kaons and pions and the influence of other dynamical contributions. 
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FIG. 1: The distribution of K/n ratio from Monte Carlo simulation for the independent production of kaons and pions with 
the global K/n ratio R g = 0.19 and mean pion multiplicity (tt) = 200. Both the distributions for kaons and pions are assumed 
to be Poissonian. The solid curve is a Gaussian fit to the distribution. 
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FIG. 2: The mean K/tt ratio from an event-by-event analysis relative to R g as a function of (n) for the cases with the Poisson 
and the MNB distributions for pion multiplicity. 
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FIG. 3: The relative variance of the distribution of K/n ratio as a function of (tt) for given R g = 0.19 for the cases with the 
Poisson and the MNB distributions of the multiplicities of kaons and pions. The experimental point from jl8| is given by a 
solid square. 



